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Abstract 

Fermion production in an external Coulomb field on de Sitter expanding universe is studied. 
The amplitude and probability of pair production in an external Coulomb field are computed 
and the cases of large/small values of the expansion factor comparatively with the particle mass 
are studied. We obtain from our calculations that the modulus of the momentum is no longer 
a conserved quantity. We find that in the de Sitter space there are probabilities for production 
processes where the helicity is no longer conserved. For a vanishing expansion factor we recover 



^^\ the Minkowski limit where the amplitude of this process vanishes. The rate of pair production in 



Coulomb field is found to be important in the early universe when the expansion factor was large 
comparatively with the particle mass. 
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I. INTRODUCTION 



The problem of pair production in an external gravitational field was intensively studied 
in the literature beginning with the Schrodinger paper [l|, who first suggested that in an 
expanding universe particle production could arise. The next important step in this di- 
rection was done by L.Parker JTHOl] who understood quantitatively the phenomenon of pair 
production in a gravitational field. The main results of L.Parker's work are related to the 
computation of the number density of scalar particles and Dirac particles created in an ex- 
panding background. Also, the limit of zero mass and infinite mass was addressed and in 
addition, an estimation of the amount of particle production for the present day expansion 
was made. Therefore in the present paper we will compare our results with these impor- 
tant predictions made by L.Parker in [7H9(], even if we will use another approach based on 
perturbations. 

In this paper we will follow the usual quantum field theory approach, were the transition 
amplitudes are defined as in [l2| using perturbation methods. This approach was also used 
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-llTJ. 20|-|22|. where transitions that generate particle production from vacuum were 
also studied on Friedmann-Robertson- Walker metrics. With this general approach, we can 
study the limit of these transition amplitudes calculated on de Sitter QED, for a vanishing 
expansion factor. This corresponds to the flat space limit, where the amplitude of this 
process vanishes and the connection with the Minkowski QED can be done. Also, from our 
general result, one could address the interesting case of large expansion, which corresponds 
to the early universe where the rate of pair production becomes important. As it was well 
established in the literature, particle production could arise as the result of the fact that the 
in and out vacua are not the same [l8|, [l9|, |27|, or because of the loss of the translational 
invariance with respect to time [12[ which in loose terms refers to the well known fact that 



in this geometry the laws of conservation for energy and momentum are lostl7H9|, ll2|. It is 

fifinn 

well known that in the Minkowski theory [5|, |6|, 113|, \1^ the pair production in a Coulomb 
field could not arise as a perturbative phenomenon. Contrary to this, in de Sitter space the 
translational invariance with respect to time is lost and one should expect a nonvanishing 
probability of transition for the process: vacuum -^ e^ + e^ in the presence of a Coulomb 
potential. It is also important to mention that the electron-positron pair could also be 
absorbed by the Coulomb field in the de Sitter geometry, which represents the time-reversed 



process. The transition amplitudes will be calculated here using only the mode expansion 
in momentum basis, which have a defined momentum and helicity. The fermion production 
in a Coulomb field in de Sitter space was to the best of our knowledge not studied so far. 
We also specify that our results are in good agrement with the work previously done on 
the particle production in a pure gravitational field, which predicts important amounts of 
particle creation in the early stage of the universe 7H9|. 



It is worth to mention that important results were obtained in 18|, ll9|], where the pro- 
duction of fermionic particles and scalar particles in electric fields on de Sitter geometry was 
considered, using a nonperturbative approach. With this approach one can obtain the den- 



sity number of particles created in the 
related to the mass of the particle 19 , 



ar out region and because there are some constrains 



271], the interesting case of large expansion from the 



early universe is not properly addressed. For the Dirac field, the results presented in 
show that there are nonvanishing rates for fermion production of mass zero and of finite 



mass in electric fields in de Sitter space. A 
to the production of zero mass fermions. In 



"ew comments need to be made here related 
7H9| and more recently in [27| it was proven 



that because in the massless case the field is conformally coupled, there is no production of 



fermions with zero mass. It will be interesting to compare our result which was obtainec 
ing the exact solutions of the field equation with the above mentioned results [7H9 



0J9, 



in the massless case. This is because the results from |7|- 
with 18 



us- 

3, 



27|, seems to be in contradiction 
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in the null mass case. 

From another point of view, it is important to directly prove that the transition ampli- 
tudes of pair creation in an external field contain terms that break the conservation laws, 
which seems to have been known for a long time but was never proven in a direct calcula- 
tion using the solutions of free field equations. For these reasons, we hope that the results 
presented here will be of interest for cosmology and quantum field theory on Friedmann- 
Robertson- Walker backgrounds. 

The second section of our paper begins with a short reminder of the work that was done 
in the theory of free fields on de Sitter QED. In the third section we present the main steps 
for calculating the transition amplitude corresponding to the process of pair production 
in the Coulomb potential in de Sitter geometry and we compute the probability for pair 
production. Physical consequences of our calculations are also analyzed in this section. The 
fourth section is dedicated to the study of limit cases as the expansion factor is small/large 



comparatively with the mass of the particle. Our conclusions are summarized in section five. 

II. FREE FIELDS ON DE SITTER SPACE 

The definition of transition amplitudes is expressed with the help of the exact solutions of 
the free Dirac and Maxwell equations on the de Sitter geometry, expanded in the momentum 
basis. For this reason, we make a short review of the results from free field theory in this 
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background. The line element [2[ for de Sitter universe is: 

ds' = dt" - e^^'dx^ (1) 

where u is the expansion factor and a; > 0. Now we know that for defining fields with half 
integer spin on curved spacetime it is required to use the tetrad fields |2[ ep{x) and e^(x), 
fixing the local frames and corresponding coframes which are labelled by the local indices 
/i, z/, ... = 0, 1,2,3. The form of the line element allows us to choose the simple Cartesian 
gauge with the nonvanishing tetrad components 

e2 = e-"*;ei = (5ie-"* (2) 

so that ep: = e'^e^ and have the orthonormalization properties e^ep = 77^, 
e^e-p = 6^ with respect to the Minkowski metric rj = diag{l^l, —1,-1). 

Now let us introduce unit normalized helicity spinors 6| for an arbitrary momentum 
vector p and denote them ^x{p) and rj^j^p) = i<J2[C,a{p)]* 

apUp) = 2pXUp) (3) 

with A = ±1/2 and where a are the Pauli matrices and p =\ p\. The particle spinors have 
the form 

P3 + P ^ \ ^ / -N P3 + P 



Then, the positive/negative frequency modes of momentum p and helicity A that were 
derived in 3|, assuming that gamma matrices in Dirac representation, are: 

U it x) = V^^ ( i^"''^- (Se--*)a(l?) , ,,_,., 
V^At X) - V^^ ( ->^--^'^'Hl^^(^^--')v.ipl , ^_....... 



where Hi, (z), Hi, \z) are Hankel functions of the first and second kind and fc = — , z/± = 
i ± ik. It is also worth to mention that these solutions for the Dirac field in the Robertson- 
Walker metric were obtained earlier in 4|, with the specification that the normalization 
constants were not determined. As was shown in [Sj, these constants depend on physical 
quantities such as momentum and play a central role in an exact calculation. So we will use 
in our further considerations the solutions from [3|. 
These solutions are normalized such that |3|: 

with the specification that the V solutions obey similar relations and the two modes are or- 
thogonal to each other. For a better explanation of what mean positive/negative frequencies 
we can use the expansion of the Hankel functions: 

For t — > — oo, the modes defined above behave as a positive/negative frequency modes with 
respect to conformal time tc = —e~^^/u: 

Up^At,^) ~ e'^^'-, V>,A(t, X ) ~ e'P'% (8) 

which is precisely the behavior in the infinite past for positive/negative frequency modes 



which defines the Bunch-Davies vacuum 12|. 



As in Minkowski space, the negative frequency modes are obtained with the charge con- 
jugation operation J5|, ll3| 

Up^Ax) ^ Vp^^x) = ^7V(f/p,A(x))^, (9) 

because the charge conjugation symmetry remains valid in a curved background. 

In our calculation, we need the form of the Coulomb potential on de Sitter spacetime 
which is dependent on the line element. The situation here is simple if we recall the conformal 
invariance of Maxwell's equations [285. The Coulomb field in Minkowski space is A^ = ^. 
Then one finds for the corresponding de Sitter potential the following formula: 

ylO(x) = ^e-"*,A^(x) = 0, (10) 

where the hatted indices indicate that we refer to the components in the local Minkowski 
frames A^ = e'^A^ . We also observe that (9) is just the expression from fiat space with 
spatial distances dilated/contracted by a factor e^"^*. 



III. AMPLITUDE CALCULATION 



In this section we will calculate the transition amplitude for pair production in a Coulomb 
field and we will explore the physical consequences of our result. The form of the transition 
amplitude can be established using the same methods as in the fiat space case, as was shown 
in 02, 
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2l| . For pair production in an external electromagnetic field, supposing that 



the fields are coupled by the elementary electric charge e, the expression of the transition 
amplitude is: 

A-e+ = -le j d^x[-gix)f'Uj,,x{xh,A^{x)Vj,',x'{x) (11) 

where e is the unit charge of the field and the fields U^,x{x) ,V^f^x'{x) are supposed to be 
exact solutions of the free Dirac equation in the momentum basis. The amplitude defined 
above for particle production in the potential fITU]) . can now be expressed with the help of 
the solutions of free Dirac equation i^. The integral that results splits in a spatial integral 
and a temporal integral. The form of the spatial integral is better known from fiat space 



theory 



a] and this result is incorporated in the next formula for the amplitude: 
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sgn{X') / dzzHl'J{pz)Hi']{p'z) 
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+sgn{X) / dzzHl^]{pz)Hl^^{p'z) 
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^t{p)v\'{p') 



(12) 



where the new variable of integration is z = —tc and sgn is the signum function. The 
temporal integral which contains the infiuence of the gravity in this process is not so simple. 
We shall give here only the main steps for solving the integrals. Using the relations between 



Hankel and Bessel functions J 
expressed in the form: 
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11| . all the integrals from the above amplitude can be 



dzzJ^{pz)Jy{p'z), 



(13) 



which are known as Weber-Schafheitlin integrals 10|, l26|]. Some aspects related to these 
integrals are given in the Appendix. 

After we complete the calculations, the final result for the transition amplitude is ex- 
pressed in terms of Gauss hypergeometric functions 2-^1 , Beta Euler functions B and unit 
step functions 9 as follows: 



j\p_ 



167r|p + p'p 



p 



-sgn{\') (p-'eip-p')u ( jj +p'-'e{p'-p)h (^^ 



p 



+sgniX) (p-'eip-p')f, ( ^j +p'-^eip'-p)U [j, 



^t{p)rix'{p').{U) 
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We specify that in the integrals (fT2|) for the case p' > p {in the case p' < p the analysis is 
similar) we let the momentum p ' have a small imaginary part p' -^ p' + ie which assures 
the convergence of our integral and makes the unit step function and the fk\^] ■, fk\ — ] 
functions well defined . In equation flT^ the newly introduced functions /fc ( ^ ) , or /^ f — j ^ 
are defined as: 

^'[p'J '^ ch^ink) B{-ll + tk) 

Uy^'\FJll + ik-l + zk-^fry-z0\ Uy\FJll-^k■l-^k■^fry~^o] 
^ ch^ink) B {-tk, I + tk) ch^lnk) B [ik, i - tk) 

ch^ink) B{l,l-tk) ' ^ ^ 

with the observation that f-k ( ^ 1 is obtained when one makes the substitution k —^ —k 
and fk\ — ] is obtained when one makes the substitution p ^ p'. To simplify our formulas, 
we can introduce the following notation x = p/p'{= p'/p) ■ A notable observation is that 
the result can be expressed without letting the momentum have a small imaginary part, 
in which case the argument x i^i the hypergeometric functions must be considered in the 
interval < x < 1- This is the domain of convergence of the hypergeometric functions, 
because in the limit x ~^ 1 the hypergeometric functions become divergent. In this case, 
one should naturally ask what happens with the result of our integrals with Bessel functions 
in the limit X = 1. In this limit, the result of the integrals can be expressed in terms of delta 
Dirac functions S{p — p'), however, when summed, these terms cancel out between them, for 
each of the two integrals with Hankel functions. 

From the final result given in flT^ and flTSj) . it is obvious that the transition amplitude 
is nonvanishing only for p ^ p' and from here we conclude that the law of conservation 
for the modulus of the momentum is lost in de Sitter space-time. One can recall that in 
this geometry the translational invariance with respect to time is lost and from here it is 
obvious that the laws of conservation for energy and momentum are lost. This means more 
precisely that the absence of Poincare invariance allows energy and momentum to appear 
from vacuum 12|. Our result related to the loss of the momentum conservation law, can 



also be seen as a consequence of the redshifting of the physical momentum in this geometry. 



Pa = e^P^ 



Pph = Pc/a(t) =pce '^ , 



(16) 



where Pc represents a comoving momentum. 

It is also more intuitive to imagine the following situation. We know that the energy and 
momenta are conserved for a particle which interacts with a constant external field. So in 
our case it is not surprising that the momentum conservation is broken due to the fact that 
the perturbation potential acts as an external time dependent source. Thereby, we conclude 
that the expansion of the spacetime will favor the production of electrons and positrons from 
vacuum with different modulus of momenta p' ^ p, in an external Coulomb field. 

It is clear that the contribution which is specific to the de Sitter geometry in our ampli- 
tude is contained in the function fkix)- This quantity encodes via the dependence of the 
parameter k = m/oj, the effect of the expansion of space on the pair production process 
and, it is the quantity that enters in our probability and help us to take a closer look at the 
limit cases. For studying the validity of our result, in the sense that the functions defined 
in flT^ are not divergent, we plot the real part (Re[f]) and imaginary part (Ini[f]) of fk{x) 
as function of the parameter k for different values of x = p/p' {p' /p)'- 
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FIG. 1: The real part of fk{x) ^^ & function of k. The solid line is for x = 0-1 ^-iid the dashed line 
for X = 0-9. 
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FIG. 2: The imaginary part of fkix) &s a function of k. The sohd hne is for % = 0.1 and the dashed 
hne for x = 0.9. 
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FIG. 3: The real part of fk{x) as a function of k. The sohd hne is for x = 0.001 and the dashed 
hne for x = 0.000001. 
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FIG. 4: The imaginary part of fkix) ^s a function of k. The solid line is for x = 0.001 and the 
dashed line for x = 0.000001. 

Our graphs show that both the real part and the imaginary part of the function fk{x) are 
finite in origin and are very convergent for large values of k (see Figs.([T])-(jl])). As x takes 
small values, we observe that these functions become oscillatory (see Figs. ([1])- ([2])). This 
oscillatory behavior is the result of the behavior of Gauss hypergeometric functions as their 
algebraic argument x^ approaches zero, combined with the oscillatory factors x^*'^- 

Because we use here the methods based on perturbations, the outcome of our calculations 
must be the probabilities of transitions. Our results for the amplitude are used in defining 
the probabilities of fermion production in de Sitter space. The main quantity that enters in 
the probability is |^e-e+| ^; so we will focus on studying its properties. We can write down 
the expression of the probability of transition by squaring the amplitude: 



V, 



^El^e 



AA' 



^E 



4/72 



e^Z 



AA 



-sgn{X)sgn{X') { fk 



^Mip)vy{p')np'-'e{p'-p) 



/-.(J^)+/.*(f;)/- 



+ 



f- 



- sgn{X)sgn{X') { fk 



fl 



fk 

p_ 
p' 

' Ik 



P_ 
pi 



/-A 



p_ 
p/ 

p'^0{p — p' 



fk 



f-k 



(17) 
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The total probability of electron-positron pair production in a Coulomb field will be: 



,tot I ^ d^P d^P' 



and via the squared amplitude |^e-e+|^ it is also determined by the factor k = m/u. 
Carrying out the integration over p ,p' is a complicated task. The form of the integrals 
that need to be solved are not known in the literature and only a numerical estimation 
can be made. However, some physical consequences can be obtained from ( 1T71) by drawing 
the graph of probability of transition Ve-e+ as a function of k, for different values of the 
ratio X = p/p' (p'/p)- As we know from Minkowski QED, the helicity conservation is a 
property which becomes important in the ultra relativistic regime. Let us see what happens 
in the not so simple de Sitter case. From the above formula for probability of transition, we 
observe an interesting property, namely that in de Sitter geometry there are nonvanishing 
probabilities for both production processes which conserve or do not conserve the helicity. 
In the process under study it is clear that helicity is conserved when A = —A', so that the 
initial zero helicity equals the sum of the final helicities A + A' = 0. Contrary, the process 
in which helicity is not conserved is obtained when A = A '. It is interesting to obtain from 
here some quantitative results from the probability formula flT7|) . These results will help 
us understand if the expansion of the spacetime will favor pair production processes with 
conserved/non-conserved helicity. 

First we represent our probability of transition fITTj) as a function oi k = m/u, for x = 0.1 
and X = 0.4, both corresponding to close values of the modulus for the two momenta p' ^ p, 
but not equal. For x = 0.9, our graphical analysis shows that the probability for a helicity 
non-conserving process is zero, and for this reason we begin our graphical representation 
with X = 0.4. 
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FIG. 5: Vf,-f,+ as a function of k for x = 0-4- The dashed hne represents the case of hehcity 
conservation and the solid line the case when helicity is not conserved. 
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FIG. 6: V (.-(.+ as a function of k for x = O-l- The dashed line represents the case of helicity 
conservation and the solid line represents the case when helicity is not conserved. 

Our numerical and graphical analysis above show that the probability for pair production 
with opposite helicities A = — A' (helicity conservation) is sensibly bigger than the probability 
for production of pair with equal helicities A = A' (helicity non-conserving case) in the 
case u > m, (see Figs. (I5])-(J6])). For m = 0, the probability of pair production in the 
helicity non-conserving case is zero, while the probability for a conserving helicity process 
is finite. We note that the zero probability of production of fermions with zero mass in the 
helicity non- conserving case is consistent with the identically null amplitude in the case when 
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helicity is not conserved due to conformal invariance. Also, the graphs (see Figs. ([5])- ([6])), 
show that when m/u increases the probabihty of pair production with opposite hehcities 
decreases and the probabihty of production of fermions with equal helicities becomes finite 
. Summarizing our result for a; > m, if the two momenta are close as moduli p ' ~ p but not 
equal, the expansion of the space will favor production processes which preserve the helicity 



conservation. This result is in good agrement with [23|] , where the helicity conservation was 
discussed in the Coulomb scattering on de Sitter space. The result from 23j proves that 
helicity conservation in Coulomb scattering is preserved as long as we deal with equal initial 
and final momenta pi = P/- 

We also see from Figs. ([SD-(IS])j that our probabilities in both conserving/non-conserving 
helicity cases, become negligible for large m/u, and from here we conclude that the produc- 
tion rate was important only in the inflationary regime, when this ratio was very small . 
Looking at the helicity non-conserving case (A = A'), we obtain from our graphical analysis 
that for X = 0.4, or equivalently as x is closer to unity, the probability for a process where 
the helicity is not conserved becomes very small. The conclusion is that in de Sitter space 
helicity conservation seems to be a property which is preserved as long as we deal with pair 
production that have the ratio of the momenta p/p ' closer to unity. 

Let us consider now the case when x is very small. This case corresponds to the situation 
in which, for example, p' ^ p (for p' '^ p the analysis is similar), so that the momentum of 
the particle p is much smaller than the momentum of antiparticle p'. We must specify that 
large momenta doesn't refer here to a relativistic momentum, but rather to the situation 
when one of the momenta is almost zero and the other has a value much smaller compara- 
tively with the relativistic case, i.e p G (10^^, 10^^) and p' G (1, 5). Plotting our probability 
flT7|) as a function of m/u we obtain the results shown in Figs. ([7])- ([9]): 
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FIG. 7: 7^6-6+ ^-s a function of k for x = 0.01. The dashed hne represents the case of hehcity 
conservation and the sohd hne represents the case when hehcity is not conserved. 
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FIG. 8: 'Pf,-e+ as a function of k for x = 0.0001. The dashed hne represents the case of hehcity 
conservation and the sohd hne represents the case when hehcity is not conserved . 
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FIG. 9: 'Pe-e+ ^-s a function of k for x = 0.0000001. The dashed line represents the case of hehcity 
conservation and the sohd hne represents the case when hehcity is not conserved . 

The first observation in this case is that the probabihty for a production process were 
hehcity is not conserved, becomes important as x ^ l? (Fig-®)- We see that in this hmit 
the pair production processes that break the hehcity conservation law become important 
as the two probabihties that correspond to the hehcity conservation/non-conservation pro- 
cesses become approximatively equal (Figj9]), for u > m. Another interesting result is that 
the probabilities for pair production in both conserving/non-conserving hehcity cases have a 
series of maxima (resonances) and minima (Fig.dT])-®), with the observation that the am- 
plitudes of these resonances decrease as the ratio m/u increases. This oscillatory behavior 
of the probability becomes more pronounced as x approaches zero. The limit x ~^ could 
help us understand the case in which one of the produced particles has a small momentum 
close to zero comparatively with momentum of the other particle. Let us make a more 
detailed analysis of this situation. It is known that in an expanding de Sitter universe the 
following first or derpro cesses could arise as the result of the loss of energy and momentum 



conservation laws 
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21| : the electron or positron can absorb or can emit photons 



(e^ ^ e^ + 7), one photon electron-positron pair creation or annihilation of the pair in a 
single photon (7^6" + e"^) and that electron-positron and photon triplet could be pro- 
duced or absorbed by the de Sitter vacuum, all effects arising as perturbative phenomena. 
The case x *^ 1 could be now explained in terms of the above mentioned reactions. One of 
the Dirac particles arising from vacuum in a Coulomb field could absorb one photon, and 
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if we take in consideration that this photon can be produced from vacuum in the above 
mentioned reactions, we can resume that the Dirac particle absorbs energy and momenta 
from vacuum. Opposite to this situation, the emission of one photon is responsible for the 
small momentum of the other Dirac particle. The emitted photon could end up in another 
electron-positron pair which could annihilate in vacuum, resulting that the momentum of 
the particle can be absorbed by the vacuum. Our remarks also help us to understand why 
the helicity is no longer conserved in the case x ^ 1- If we start, for example, with an elec- 
tron, it could emit or absorb one photon and the spin in this reaction is no longer conserved. 
As a final conclusion, in the limit of small x ^6 establish that helicity conservation is no 
longer a rule and the processes where helicity is not conserved become equal in importance. 

It is also important to mention that our graphical results for probability (Figs.([S])-([n])), are 
valid for small momenta (p, p'). This is because the factor \'p + p'\^^ cancel our probability 
in the limit of large momenta. Thus we can conclude that particles with low momenta were 
produced in external Coulomb field, due to the early expansion of the space. This is in 
agreement with the formula for physical momentum (IT6ll . from which we see that for large 
expansion factor uo the physical momentum become small. 

As a general observation from all the graphs of probability as a function of A; = m/oj 
it is clear that the amount of particle production was important only for small k. This 
corresponds to the early universe, when the expansion was large comparatively with particle 
mass (w ^ m). This result is in good agrement with the previously obtained results [7- 
ISj, which also predicts that the amount of particle production in the early universe was 
important. 

Taking in consideration that the particle production rate was important for small values 
of fc ^ 1, it will be interesting to discuss this case. This is the case of a large expansion 
factor, with non-zero mass. The amplitude and probability formulae in this section are used 
in our analysis and, in addition, we have at our disposal the graphs for probabilities. We 
specify that the two cases for conserving/non-conserving helicity are considered. For our 
considerations it will be sufficient to restrict ourselves, for example, to the case p < p'. We 
use helicity bispinors, for A = A ' and A = —A ', and restrict the analysis only to the square of 
our amplitude. In an orthogonal local frame {cj} we take the electron and positron momenta 
in the plane (1,2), denoting their spherical coordinates as p= {p,a,(3) and p' = {p','y,{p), 
where a, 7 G (0, vr); /3, </) G (0, 27r) . Then ioi (3 = tc ,ip = and p < p' the final result for 
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the square amplitude is: 



e^7? 1 

\A ,|2 = 1^ ^_Q(j^> _y\ 



h\j, 



2 



f-k \ — 
pi 



+ 

cos2 (^) for A = -A' 



Sin2 /a±7 



(19) 
^) for A = A' 



From this formula, if we set a = vr , 7 = 0, corresponding to the case when the momenta 
of the electron and positron are parallel but move in opposite directions, we obtain zero 
probability for a helicity conservation process, while the probability for a process where 
helicity is not conserved, becomes maximum. All the particular cases where the momenta 
are oriented on different directions can now be obtained from fll9p . We can conclude that for 
A; ^ 1, the expansion of the space will favor helicity non-conserving processes in which the 
electrons and positrons have parallel momenta and move in the opposite directions. This 
result shows that only in the helicity non- conserving processes the chances of separation 
between electrons and positrons in the early universe become important. Thus it seems that 
for a = 0,7 = the situation is opposite, conserving processes being dominant, with the 
specification that this case corresponds to the pair having parallel momenta and moving in 
the same direction, which increases the probability of annihilation of the pair. 

In the next section we shall pay a closer attention to the analysis of this probability 
in the limit cases of the parameter k. It will be interesting to study the probability in 
the limit of large k to make a connection with the well established results from literature, 
where the transition from vacuum to non-vacuum states are studied using the Bogoliubov 
transformations. This is because we expect that the behavior of our amplitude/probability 
for large k to somehow reproduce up to some factors, the results obtained when one uses 
the asymptotic solutions of the field equations. 

IV. LIMIT CASES 

The de Sitter metric reduces to the Minkowski one when the expansion factor of the space 
vanishes, w = 0. In this limit the free field equations (Dirac and Maxwell equations) also 
become the field equations from fiat space, so we expect our amplitude to reduce in this 
limit to the corresponding amplitude from Minkowski space, which is zero. In the fiat space 
this process of pair creation from Coulomb field is forbidden by the laws of conservation for 
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energy and momentum. The Minkowski limit corresponds in our case to /c = oo, which is 
the ratio (mass of particle/expansion factor). 

Let us see first what happens for large values of the parameter k. In this case one can 
observe from fllSp that the second and the third arguments in the hypergeometric functions 
2-Fi(a, b; c; z) become approximatively equal (6 ~ c) and we can use 2-^1(0, b;b; z) = (1 — z)~°- 
to simplify these functions. For Beta Euler functions one can apply the Stirling formula to 
obtain an approximation in the limit of large k: T{z) = z^~^^'^e~^\/27i{l + 0{z~^)), \z\ — ?■ cxd. 
Finally, one finds for large k (only the terms proportional with e'^^ are considered in the 
functions fk, f-k)- 

-. {l-X^yl- '^-^^) - V^ 2(l-x^)3/2 ' (20) 

with the specification that the term e^'^^ originates from ch^'^{nk). We see that these 
functions are highly convergent and vanish as: 

f±k{x) ~ e--^ (21) 

"or k — >■ 00. An important observation is that the factor e"'^'^ is used in the literature 
a y y a ^^ cWcten.„g t^e vacuu. .o no„-vacuu. .an.t.n. i„ .h. geon,e... 
As noted before, the well established results from the literature use the asymptotic behavior 
of Hankel functions to establish the positive and negative modes and then with the help of 
the Bogoliubov coefficient Pp, one can obtain the number density of particles per volume 
|/3ppT2^. For the above mentioned vacuum non- vacuum transitions [12, llSl, ll9|, |27|, the 
Bogoliubov coefficient is proportional to |/3pp ~ g~27rm/a;_ This exponential factor is also 
found in our probability for large k. Using equation (120|) one can obtain the probability of 
pair production in a Coulomb field for k ^ 1 (or m ^ u). The cases p < p' and p > p' 
are considered in our analysis, with the observation that we no longer need to use step 
functions. Using the spherical coordinate oip= (p, a, vr) andp ' = (p',7,0), the final result 
for the probability (conserving and non-conserving helicity processes are included), in the 
case m ^ oj, is: 

p ^ e'Z^ e-'-^ PP' [ cos^ (^) , for A = -A' 

''^ 2^Q7:^\p + pr{p"-p^f\sin^{^),i0T\ = \'. 



Recently, in 27|, it was shown that the condition m ^ u must be imposed in order to 



obtain production of a thermal spectrum of radiation, as calculated in the work done by 
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Gibbons and Hawking [25|. The exponential in the probabihty equation (122|) has the form of 
a Boltzmann factor, that means that the spectrum of particles created in our case for m ^ u 
is thermal, having the Gibbons-Hawking temperature T = uj/27t. Also^t is remarkable that 



our probability is proportional with e"^'^™/'^, which was found in 12|, |l8|, [l9|, |27|, to be the 
factor which characterizes the probability of particle production from vacuum for m ^ u . 
Now let us make an estimation of this probability for the present day expansion of the 
Universe. The last astronomical observations give u = 2,3 ■ lO^^^s"^, for the Hubble 
constant. Then from the ratio mc^/hw ~ 5 ■ 10^"^, (where m is the electron mass) one can see 
that the large mass approximation is good in considering the pair production at present. It is 
immediate from (122]) . in which we replace the above mentioned ratio, that the probability of 
pair production at present is very small, and we can say that there is no fermion production 
of finite mass in Coulomb field due to present day expansion, a result which is in good 



Qlsi. 



agrement with the conclusion from 

We can now address the problem of the Minkowski limit. The graphs (Fig. ([T])- (14])) from 
the previous section, for functions fk{x)i invite some comments. First, we see that in 
the limit /c — ;■ oo, the functions fk{x) vanish, which is in agrement with the result ob- 
tained analytically, fooix) = 0- From here we recover the Minkowski limit, where particle 
production in a Coulomb field is forbidden by the laws of conservation for energy and 
momentum. As k goes to large values the de Sitter amplitude equal the Minkowski one 
Ae-e+\k=oo = ^e-e+(/^«'^) = 0, which corrcspouds to the fact that in an asymptotic flat 
universe these amplitudes vanish. Also it is worth mentioning that for large momenta, 
p,p' — 7- oo, the amphtude (JH]) and probability ( TT7]) vanishes and the recovery of Minkowski 
results would probably not be surprising. 

Secondly, it is obvious from our graphs (Fig.(IT])-(l4])) that our formula for functions fk{x) 
becomes important when the expansion factor is larger than the mass of the particle. The 
null mass limit corresponds to the early stage inflation of the universe, when the expansion 
factor was large comparative to the mass of the particle w — t- oo. The fk{x) function 
simplifies in this case (i.e. 2Fi{l, 1/2; 1/2; z) = 2^1(1, 3/2; 3/2; z) ~ (l-z)-^) and for A; = 
we obtain: 

fo{x) = -a+xr'. (23) 

Also, it is immediate from (JH]) that for fc = the amplitude of transition is nonvanishing 
only for A = — A'. This is in agrement with the observations made in the previous section, 
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that the most probable transitions in the null mase case are those in which the helicity is 
conserved. Replacing f l23|) in the transition amplitude flT^ and observing that in this case 
we no longer need to use unit step functions, the final result is: 

Att'^IP + p '\^ p + p' — tO 
It is remarkable that the amplitude of transition in this limit depends only on the modulus 
of two momenta p, p' . The same results are obtained from our integrals in ( fT2|) if we set 
A; = 0, in which case, the Hankel functions simplify to: 



r(2), 



2 



1/2 



and the integrals from f lT2|) can be evaluated adding a factor e""^^, where the parameter e 
is positive, to assure the convergence in the limit z — )■ oo. Note also that the amplitude in 
the zero mass limit vanishes if the particle and antiparticle have the same helicities, which 
corresponds to the non-conserving helicity in the process presented above, which also confirm 
our graphical results for probabilities. This is an expected result since the Dirac field for a 
null fermion mass is conformally invariant and the de Sitter metric is conformal with the 
Minkowski one. 

One can also see that for a very large expansion factor but non-zero mass, the exponen- 
tial factor e~'^*, from the Coulomb potential, implies that only negative large time will be 
relevant, which means that we can consider z ~ e~^^ large everywhere in the integrals ( TT2l) . 
Translated in our calculations, this implies that we can use the expansion: 

/ 2 \ ^/^ 
Hf\z) =[ — ) e-'(^-'^"/2-V4)^ ^^^ (26) 

which differs from H\^L{z) only by the phase factor e«(M'r/2+7r/4) _ ^\^^q conclusion is that for 
very large expansion factors, the amplitude of this process will be mainly determined by the 
contributions from the distant past when space was contracted. It is from here immediate 
that the frequency blueshift in the distant past makes the mass irrelevant, a fact which leads 
back to a Minkowski conformal theory. 

Further, let us consider a closer look at our amplitude in the null mass case ^IM . As 
before, kinetic parameters can be introduced in the orthogonal local frame {cj} . We take 
the particle and antiparticle momenta in the plane (1, 2) denoting their spherical coordinates 
by P = (p, tt,/3) and p' = (p',7,(/?) where a, 7 G (0, tt) and (3, ip E (0,27r). Then, using 
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the helicity bispinors, one can obtain for A = — A' = | the final expression of the transition 
amphtude: 

e^Z 1 



j\p 



l(fe=o) 



47r^|p + p 'P p + p' — iO 



coMf)cos(I)+e*'-''>sm(|)sm(2 



(27) 

where the terms in the square brackets were obtained from C,^/2{p)V-i/2{p ')■ From the above 
formula, one can obtain all the particular cases of interest in which the electron and positron 
move in the same direction or move in opposite directions. The probability of production of 
null mass fermions, if we set /3 = , y? = vr, will be: 

From the above formula we see that for a = vr , 7 = our probability vanishes. This 
corresponds to the case in which the particle and antiparticle have parallel momenta and 
move in opposite directions. Now if we set a = ,7 = the probability takes it's maximum 
value, and this is the case in which the particle and antiparticle momenta are parallel and 
they move in the same direction, increasing the probability of annihilation of the pair. The 
obvious conclusion is that if the null mass fermions are produced they will rather annihilate 
each other, than become separate. From here we conclude that there are small chances of 
production of null mass fermions in the early universe. As a final remark, one can see that 
our perturbational calculations in the null mass case is closely to those from the literature 
SU, wMch e.ab„.K. that t.e. . no p.oduc^on of .«o ... fe„.on. ,n t.e e„l. 
universe. 

V. CONCLUSIONS 

Fermion pair production in a Coulomb field in a de Sitter expanding universe was consid- 
ered in our paper. We found the expression of the differential probability for the transition 
vacuum — )■ e^e^ in the potential A^. We considered the initial and final states of the field 
as exact solutions of the free Dirac equation in de Sitter space with a defined momentum 
and helicity. We also found that the amplitude of pair production depends, in an essential 
way, on the parameter k = m/u. For a vanishing k we recover the expected result due to 
the conformal invariance of the theory in the massless case. The final result proves that the 
law of conservation for the momentum is lost in de Sitter space as a result of the loss of 
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translational invariance with respect to time. From our amplitude formula, we recover the 
Minkowski limit where there is no particle production in the Coulomb field. Our final results 
related to helicities show that the helicity conservation processes seem to be dominant in 
the case when the momenta of the electron and positron are close in modulus (x — ?■ 1)- 
When the ratio of the two momenta is small we found that the processes where helicity is 
not conserved become important. Our study also led to an important result related to the 
mechanism of separation between matter and antimatter in the early universe, namely that 
only in the helicity non- conserving case there seems that are nonvanishing probabilities for 
producing pairs which move in opposite directions, increasing in this way the possibility of 
separation between particles and antiparticles. 

A notable result is that only particles with small momenta were produced in Coulomb 
field in the early universe. 

Another result is that for the present day expansion, the amount of fermion production in 
Coulomb field is negligible. Contrary to this, our calculations and graphs shows that in the 
early stage of the universe the amount of pair production was important. This confirms the 
results from literature related to the phenomenon of pair production in the early universe. 

Our study doesn't take into account the problem of wave packets in the de Sitter geometry. 
This problem could be important in studying the effects of pair production, and receives 
little attention in the literature. However, the study from [24] is a first step in this direction, 
where the construction and propagation of gaussian wave packets on de Sitter geometry was 
discussed both for scalar and Dirac fields in the null mass case. The study of wave packets 
in this geometry could improve our understanding of scattering processes and could help us 
to understand the physical significance of the cross section. 

VI. APPENDIX 

The main steps leading to our amplitude JHM will be detailed in this section. First, we 



will use the well known formula for Hankel functions in terms of Bessel functions J 



10|,lll|: 



Hi^\^) 



ism{nfi) 



H^'\z) = '''^'^^^'\ ^-^^'\ (29) 

^ i sin(7r/i) 



This transformation finally leads to integrals of the type Weber-Schafheitlin 10|, |26 |. 
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These integrals can be solved in two distinct cases: 

dzz-'J^{az)J^{f3z) - y 2 ! 



/u + z/ — s + 1 li — u — s + 1 o? 

''-'''' [ 2 ■ 2 ' " + 1' # 

/3 > a > 0, Re{s) > -1, Re{fi + z/ - s + 1) > 0, (30) 



and 



dzz-'J^{az)J^{Pz) f^ y 2 ) 



2^a'^-«+ir(i/ + i)r( ^-'-+^+i ) 

^ /yu + z/-s + l z/-/i-s + l /32~ 

2 1^ 2 ' 2 ' 'a2 ' ' 

a > /3 > 0, Re{s) > -1, i?e(/i + z/ - s + 1) > 0, (31) 

both cases being fulfilled in our analysis. 

Now one sees that our integrals in (iT3ll demand s = —1. This is problematic because 
the integral becomes oscillatory for 2; — )► oo. Then a simple way to solve this problem is to 
consider a parameter s of the form 

s = -1 + e (32) 

and let, in the end, e — > 0. No notable differences appear when evaluating the integrals in 
flT^ directly with s = —1 and with 11, v = 1/2 ± ik . 
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